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The associated production of a Higgs boson and a top-quark pair, tiH, in proton-proton collisions 
is addressed in this paper for a center of mass energy of 13 TeV at the LHC. Dileptonic hnal states 
of tiH events with two oppositely charged leptons and four jets from the decays t —>■ hW^ -A bl^ve, 
i -A bW~ -A bl~V(_ and h -A bb, are used. Signal events, generated with MadGraph5_aMC@NLO, 
are fully reconstructed by applying a kinematic fit. New angular distributions of the decay prod¬ 
ucts as well as angular asymmetries are explored in order to improve discrimination of tiH signal 
events over the dominant irreducible background contribution, tibb. Even after the full kinematic 
ht reconstruction of the events, the proposed angular distributions and asymmetries are still quite 
different in the tiH signal and the dominant background (tibb). 


I. INTRODUCTION 

On July 4th 2012 the ATLAS |T] and CMS [2] collab¬ 
orations announced the discovery of a scalar particle at 
CERN’s Large Hadron Collider (LHC). This new particle 
with a mass of about 125 CeV was later identified as the 
Higgs boson, responsible for the generation of all particle 
masses through the mechanism of spontaneous symme¬ 
try breaking [3]. So far, the measured properties of the 
Higgs boson have shown a remarkable consistency with 
those predicted by the Standard Model (SM) of particle 
physics. Nevertheless, it is by now clear that the SM can¬ 
not explain all the observed physical phenomena, as for 
instance it fails to provide a candidate for dark matter 
or a means to explain the matter anti-matter asymmetry 
in the Universe. However, as more data is being accu¬ 
mulated and analysed at the LHC, it becomes increas¬ 
ingly clear that any new physics theory has to resemble 
very much the SM at the electroweak scale. In the first 
run, the ATLAS and CMS collaborations have studied 
in great detail the main four Higgs production modes at 
the LHC |1], namely gluon fusion (including bb fusion). 
Vector Boson Fusion (VBF), associated production {VH, 
with V = W,Z) and tiff production, with centre of mass 
energies of 7 and 8 TeV. For each production mode sev¬ 
eral decay channels were considered and analysed in great 
detail. 

The production of the Higgs boson in association with 
a top quark-anti-quark pair, pp -A tiff |5], constitutes 
the only way (together with single top plus Higgs which 
has an even smaller cross section) to directly probe the 
top-quark Yukawa couplings. Moreover, it is also con¬ 


taminated by a huge background coming mainly from 
pp ^ ti -\- jets. For this particular production process 
several decay channels have been studied [SHE]- The very 
complex final states, together with the huge backgrounds 
to the process, makes it the most difficult Higgs chan¬ 
nel to study at the LHC. Nevertheless, with just a few 
events, both collaborations have reach a sensitivity down 
to about 2-3 times the SM value which constitutes a 
remarkable achievement. 

The current studies EHH], use the kinematic informa¬ 
tion of the events to separate the signal from the back¬ 
grounds. In this work we advocate the introduction of 
new variables that make use of the information from (lack 
of) spin correlations in the signal and background pro¬ 
cesses mini: the top and anti-top quarks are natural 
spin analysers of this process. We will show that part 
of the spin information that is present in the matrix el¬ 
ements survives the parton showering, detector simula¬ 
tion, event selection and event reconstruction. These new 
variables could play an important role in background dis¬ 
crimination, possibly leading to an improvement in the 
precision of the measurement of the top-quark Yukawa 
coupling. Even though we will consider only the ir¬ 
reducible tibb background, without a highly-optimized 
event-reconstruction method, we will argue that our find¬ 
ings are also valid in a more general and realistic case. 

We should also note that the type of variables pro¬ 
posed in this work can then be used to probe the CP 
nature of the top-Yukawa coupling mi. In many models 
like the CP-violating two-Higgs doublet model [T2] (the 
status of this model after the LHC run 1 was recently 
presented in m), CP violation appears explicitly in the 
Higgs sector via mixing of CP-even and CP-odd states. 
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The determination of the CP nature of the Higgs boson 
and its interactions is of the utmost importance at the 
LHC. Finally one should note that tiH production can be 
studied at future linear colliders such as the ILC, which 
will lead to a tremendous improvement in the precision 
of measurements of the Yukawa couplings [14] . 


II. SIGNAL AND BACKGROUND 
GENERATION AT THE LHC 

Given that the goal of this work is to study how well 
spin information can be used to improve the current 
search strategies, we only consider the signal and its 
dominant irreducible background. The signal (tiH) and 
background (tibb) processes were generated, at leading 
order (LO), using MadGraph5_aMC@NLO (TS] with the 
NNPDF2.3 PDF sets [TB]. The full spin correlations in¬ 
formation of the t bW~^ —?> bi^ve, i —>■ bW~ bi~Vi 
and h ^ bi) decays, with G {e^,/r^}, was retained by 
allowing MadSpin uni to perform the decay of the heavy 
particles. Although other decay modes of the top quarks 
and Higgs boson could be considered, in this paper we 
focus on the most challenging leptonic decay channel i.e., 
the dileptonic decay of the ti system together with a 
Higgs decaying to the dominant SM decay (H —>■ bb). 
We argue the full kinematical reconstruction of the un¬ 
detected neutrinos in such events, even if difficult, still 
preserves the angular distributions that could help in 
discriminating signal from irreducible backgrounds. The 
events were generated for the LHC with a centre of mass 
energy of 13 TeV with the default dynamic factoriza¬ 
tion and renormalization scales, setting the masses of 
the top quark and the SM Higgs boson to 172.5 GeV 
and 125 GeV, respectively. We do not attempt to in¬ 
vestigate possible departures from the SM nature of the 
Higgs boson in this paper, assumed to be a scalar parti¬ 
cle (CP = 1). The generated events were then passed to 
PYTHIA 6 (17] for shower and hadronization. In order to 
obtain more realistic results, for example for differential 
cross sections and efficiencies, we passed the generated 
events through Delphes m to perform a fast detector 
simulation of a general-purpose collider experiment at 
the LHC. We used the ATLAS default card for the sim¬ 
ulation and performed the analysis of the generated and 
simulated events with MadAnalysis 5 m in the expert 
mode [20] . 

The full kinematical reconstruction of tiH events is 
very challenging in the dileptonic decays of the ti sys¬ 
tem, since both undetected neutrinos need to be recon¬ 
structed. In this paper, we explore the advantages of fully 
reconstructing the tiH system in the dileptonic topol¬ 
ogy, by applying a kinematic fit to the events using mass 
constraints and energy-momentum conservation. Events, 
after detector simulation, are accepted if they had at 
least four reconstructed jets and two charged leptons with 
transverse momentum px > 20 GeV and pseudo-rapidity 
ry < 2.5. No cuts are applied to the events transverse 


missing energy (Pt)- 

In the following we will normalise the distributions of 
signal and background to equal area, irrespective of their 
(fiducial) cross sections and efficiencies: for our goal it 
suffices to show which observables are sensitive to the dif¬ 
ference in spin information in the signal and background 
events and how well this information can be retained in 
a realistically reconstructed event. For the same reason 
we also abstain from performing a careful analysis of the 
uncertainties in the event generation as well as the in¬ 
clusion of next-to-leading order corrections in the strong 
coupling. 

HI. RECONSTRUCTION OF DILEPTON ttH 

EVENTS AFTER DETECTOR SIMULATION 

As previously stated we will perform tiH event recon¬ 
struction in final states with two charged leptons and at 
least four jets, after Delphes simulation. We do not at¬ 
tempt to tag the flavour of jets from the hadronization 
of b quarks, i.e., we do not use any 6-tagging tool to 
help with the identification of the heavy flavour compo¬ 
nent of jets; a task left outside the scope of this paper. 
The full kinematical reconstruction requires the knowl¬ 
edge of the jet and charged lepton momenta, together 
with the transverse missing energy. We use the W and 
the top quark masses as constraints. The Higgs boson 
mass (rriH = 125 GeV) is used to maximise the probabil¬ 
ity of the best combination of two jets chosen among the 
ones which were not utilised in the ti system kinematical 
reconstruction. The transverse missing energy is re-fitted 
to improve the resolution of the experimental measure¬ 
ment. After applying the constraints, six unknowns need 
to be fully reconstructed in the dileptonic tiH events, 
which are the 3-momenta of the two neutrinos present 
in the events. To find a kinematic solution we assume 
the neutrinos are responsible for the missing transverse 
energy, i.e., 

Pl+Px=fx, (1) 

P''y+Pl=$y. ( 2 ) 

In addition we apply the following mass constraints to 
the ti system of the events. 


{pe+ + Pl^)"^ — 

(3) 

(pi- + Pv)'^ = 

(4) 

(Pw+ +Pb)'^ = TOt, 

(5) 

(pw- PPlf = ‘ml- 

(6) 


While f!x and represent the x and y components of 
the transverse missing energy, and (pt, and pi) 
correspond to the two lepton (two 6-jets) four momenta, 
respectively, from the t (i) decays, mw and rrit are the 
W-boson and top quark masses, respectively. The mass 
of the VF-boson was set to 80.4 GeV. 

We study the performance of the reconstruction with 
respect to the generated parton-level Monte Garlo in- 
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formation. To make sure the kinematical reconstruc¬ 
tion produces sensible results, the reconstruction is first 
applied to truth-matched objects, i.e., jets and leptons 
which are matched to their parton-level generated quarks 
and charged leptons, using a AR criterion (the minimum 
distance in the pseudorapidity-azimuthal angle plane, 
AR, between the reconstructed jet or lepton and the 
parton-level quark or charged lepton, ensures the match¬ 
ing) . Even though we use a rather simple kinematical re¬ 
construction method the efficiency using truth-matched 
objects is 62%. In figure the neutrino (left) and an¬ 
tineutrino (right) pt from signal events are shown. The 
generated distributions (filled histograms) are compared 
with the truth-match reconstructed ones (solid lines). In 
the bottom plot, the ratio between the two distributions 
is shown. Although a slight slope is visible in the ratio 
plot, more significant at high pr due to radiation effects 
not explicitly corrected for at the moment, good agree¬ 
ment between the reconstructed distributions is observed 
with respect to the parton-level neutrino distributions, 
making clear that the full kinematical reconstruction of 
tiH events is possible. In figures and the pr dis¬ 
tributions of the t (t) quarks and {W~) bosons are 
shown, respectively, for the tiH events. Once again we 
see a similar behaviour as observed for the neutrino pt 
distributions, i.e., good agreement between the recon¬ 
structed kinematic distributions and the corresponding 
ones at parton level, in spite of the slight slope for higher 
Pt values, in the ratio plot. Although the kinematical 
fit can correct, to a large extend, the effects of radiation, 
at high values the differences between reconstructed and 
generated distributions may require an additional correc¬ 
tion. Even though this would not difficult to implement, 
we have decide not to apply it here once it may depend 
on the exact experimental environment conditions and 
does not contribute significantly to the main discussion 
of the paper. 

In a second step, the truth match condition is dropped, 
bringing the analysis closer to what can be done at col¬ 
lider experiments. For this particular case we perform all 
possible combinations of reconstructed jets and charged 
leptons (after detector simulation) in order to reconstruct 
the top and anti-top quarks, together with the Higgs bo¬ 
son. For the ti system reconstruction we used the same 
procedure based on equations eq. We calculate 

the probability Pjf that the event is compatible with 
the equations, using probability density functions for the 
neutrino and anti-neutrino px distributions, the top and 
anti-top quarks mass distributions as well as the IT+ and 
W~ bosons mass distributions, obtained at parton level. 
To identify the two-jet combination, among the ones not 
used in the ti reconstruction, that best matches the jets 
from the Higgs boson decay, we associate to each combi¬ 
nation, a weight Ph, 

Ph = l/l^y {Pi +P]Y - mH\, (7) 

related to how close the Higgs boson mass {rriH = 
125 GeV) is to the invariant mass of each particular jet- 


pair combination. The solution with highest Pti x Ph 
is chosen as the right one for the full kinematical recon¬ 
struction of the events. This fixes completely the as¬ 
signment of jets and charged leptons to their parent t, i 
quarks and Higgs boson. Due to the increase in the num¬ 
ber of possible combinations which can satisfy eq. 0-0. 
88 % of all events are reconstructed by the kinematic fit. 
This will obviously lead to an increase of the combinato¬ 
rial background but, as we will see later, the kinematics 
are in most cases distinct from the right combinations. 
In figures and we show the pr distributions of the 
top quarks and W bosons, respectively. The kinemat¬ 
ically reconstructed px, with no jets and leptons truth 
match, is compared with the parton-level distribution. 
We see a good correlation between the kinematically re¬ 
constructed distributions with respect to the parton-level 
ones, thus ensuring that the reconstruction works fairly 
well. We did not attempt to further optimise the event 
reconstruction because, again, the main goal here is to 
show that a reconstruction is possible with a reasonable 
efficiency. 

IV. ANGULAR DISTRIBUTIONS 

We will focus on angular distributions in fully recon¬ 
structed tiH events involving three-dimensional angles 
between the decay products of the tiH dileptonic final 
states. Following the full reconstruction of events, we 
define two reference frames: 

• Frame 1: the full tiH centre-of-mass system, built 
by using the laboratory four-momenta and, 

• Frame 2: the iH centre-of-mass system recoiling 
against the t quark, in the tiH system (i.e., in 
Frame 1 as defined above). 

For the generated distributions (with and without the 
px and 77 cuts applied in the event selection), we use the 
parton-level four-momenta of all relevant objects. For 
the reconstructed distributions (with and without truth 
match), we use the four-momenta obtained after apply¬ 
ing the kinematic fit reconstruction. We define the an¬ 
gle between the Higgs momentum direction (in the iH 
centre-of-mass) with respect to the iH direction (in the 
tiH system) as 0*^ and the angle of the Y top quarks 
or Higgs decay products (IT+,IT “,b and b jets) 
momentum (in the Higgs centre-of-mass system) with re¬ 
spect to the Higgs direction (in the iH system) as Oy ■ 
We should stress the fact that, when boosting Y to the 
centre of mass of the Higgs boson, the laboratory four- 
momenta were used (in a direct, rotation-free boost). 

In figure we show distributions at parton level, with¬ 
out any cuts, for the product of cos ( 6 *^), and cos (Oy), 
for Y = £+(left) and Y = (right). We can see the 
distributions are quite different between signal and back¬ 
ground events. The effect of applying the px and rj cuts 
to jets and leptons is seen in figure A clear reduc¬ 
tion on the number of events is observed due to the cuts 
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applied. In figure we can see the effect of the kine¬ 
matic fit reconstruction, still with the truth match infor¬ 
mation. The information on the angular distribution is 
preserved to a large extent, even after the full kinematical 
reconstruction. As we will see this is also true when the 
reconstruction is performed without truth match. In fig¬ 
ure we show the reconstructed product (without truth 
match) of cos(0)|^) and cos (Oy), for Y = f+(left) and 

we show the same 


Y = £ (right). In figures 10 


and 
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distributions but with the charged leptons replaced by 
the kk-bosons and 6-quarks from the Higgs decay, re¬ 
spectively. It is quite apparent that some of the angu¬ 
lar distributions allow discrimination between signal and 
background even after the full kinematical reconstruction 
without the truth match. Since we did not try to opti¬ 
mise the kinematical reconstruction, it is foreseeable that 
better results could be obtained in the future. 


V. FORWARD-BACKWARD ASYMMETRIES 

Based on the angular distributions introduced in the 
previous section, we propose to use several forward- 
backward asymmetries (A))^) in this paper, defined 
using the double angular product 

xy = cos X cos {Oy)- ( 8 ) 

The asymmetries can be easily calculated, both at parton 
level and after the kinematic fit reconstruction, and are 
defined as, 

y _ N(xy > 0) - N(xr < 0) 

N(xy > 0 ) + N(xy < 0 ) ’ ^ ^ 

where N{xy > 0) and N{xy < 0) correspond to the 
total number of events in the corresponding angular dis¬ 
tribution with xy above and below zero, respectively. 
These asymmetries can be quite different between the 
signal tiH and the irreducible background ttbb. In Ta¬ 
ble [T] we present the values of the asymmetries, with no 
cuts applied to the events, at parton level (and at LO) 
for different choices of the final state particle (Y) that is 
boosted to the centre of mass of the Higgs boson. As we 
can see, there are clear differences for some of the asym¬ 
metries i.e., (6 from t), between 

signal and background. We show in figure an exam¬ 
ple of two-binned angular distributions for Y = £~ and 
y = €+, respectively, evaluated at parton level without 
any py or ry cuts applied to the events. 

In Table |T] we also show the values of the asymmetries 
after all cuts and the kinematic fit reconstruction (with¬ 
out truth match), for different choices of the final state 
particle (Y) boosted to the centre of mass of the Higgs 
boson. As we can see, even after the kinematical recon¬ 
struction there are clear differences for some of the asym¬ 
metries i.e., A))g^ (6 from t), between 

signal and background. Note that the two asymmetries 
A|)g^ (6 from H) and (6 from H) are zero at the 


(Asymmetries @ LO) Parton level Reconstruction 



ttH 

ttbb 

ttH 

ttbb 

^FB 

-0.157 

-0.137 

-0.141 

-0.268 

^FB 

-t0.291 

-t0.056 

-t0.331 

-t0.118 

Ay=w+ 

^FB 

-0.154 

-0.119 

-0.119 

-0.275 

aY=w- 

^FB 

-t0.317 

-t0.067 

-to. 348 

-t0.127 

A^j-F (6 from t) 

-0.155 

-0.141 

-0.179 

-0.306 

A'fF (6 from i) 

-tO.293 

-t0.053 

-tO.334 

-t0.117 

A^F (6 from H) 

-to.ooo 

-tO.OOl 

-t0.086 

-0.048 

AfS*” (6 from H) 

-to.ooo 

-0.001 

-0.086 

-t0.048 


TABLE 1. Values for the asymmetry for ttH and tibb events 
at the LHC. The second and third column show the observed 
asymmetries at the parton level (without any cuts), while the 
fourth and last column show same asymmetries after applying 
the selection cuts and the kinematical reconstruction (without 
truth match). 


parton level, but non-zero at the reconstructed level due 
to a n on- perfect reconstruction of the event. We show in 
figurej^the two-binned angular distributions for Y = l~ 
and Y = £+, respectively, obtained after all cuts and full 
kinematical reconstruction of events. Although distor¬ 
tions (that may be corrected for) are visible as a conse¬ 
quence of the cuts applied and kinematical reconstruc¬ 
tion, some of the angular distributions and asymmetries 
show significant differences between the signal and domi¬ 
nant background, even after reconstruction (see Table |T|. 

One last comment is in order in what concerns the re¬ 
constructed mass of the Higgs boson. Even after the full 
kinematical reconstruction and possible contamination 
from the combinatorial background arising whenever the 
reconstruction is performed without truth match, it is 
still possible to recognise, in the m^i variable, the mass 
peak corresponding to the right combination of 6-quarks 
coming from the Higgs boson. In figure we show a fit 
of the Higgs mass in signal events, just to guide the eye, 
performed with RooFit m using a Chebychev poly¬ 
nomial (to parametrise the combinatorial background) 
and a Gaussian distribution (to describe the Higgs 
mass reconstructed from two 6-quarks). Once again no 
optimisation is performed in the fit. The effect of the 
combinatorial background is clearly visible as a shoulder 
towards lower values of the invariant mass distribution 
which extends to higher values with a long continuous 
tail. We argue that it is important to understand the 
different components of the combinatorial background 
and dedicated studies must be performed to minimise 
the effect of its uncertainties, but this is largely outside 
the scope of this paper. 









5 


VI. CONCLUSIONS 

In this paper the tiH production in proton-proton col¬ 
lisions at the LHC is addressed, for a centre of mass 
energy of 13 TeV. Fully reconstructed, dileptonic final 
state tiH events, from the decays t —>■ bW~^ —)■ 
i —)■ bW~ —t U~vi and h —)■ bb, are used to probe new 
angular distributions and asymmetries that allow bet¬ 
ter discrimination between the signal and the main irre¬ 
ducible background. We show that it is possible to fully 
reconstruct tiH final states in the dileptonic topology 
and, even with a reconstruction which is not optimised, 
still be sensitive to the new angular distributions and 
asymmetries, which seem to be quite different between 
the signal and background even after full reconstruction. 

One should again stress that current experimental re¬ 
sults on the pp —>■ tiH are already very impressive even 
though, essentially, kinematic variables are used. We 
have shown that the use of new variables that make use 
of the spin information of signal and background pro¬ 
cesses can further improve the results for the cross section 
measurement. Furthermore, the spin information that is 


present in the matrix elements survives showering, detec¬ 
tor simulation, selection and reconstruction, even in the 
most challenging decay channel of dileptonic tiH events. 
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FIG. 1. Neutrino (left) and antineutrino (right) pr distributions. The generated distribution (shadowed region) is compared 
with the kinematical £t reconstruction with truth match (full line) distribution. 
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FIG. 2. Same as in Hgure[^ but for the top (left) and anti-top quarks (right) pr distributions. 
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FIG. 3. Same as in figure. but for the (left) and W (right) pr distributions. 
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FIG. 4. Reconstructed top (left) and anti-top (right) quark pt using the kinematical fit (without truth match) as a function of 
the Pt at parton level. 
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FIG. 5. Same as in figure. but for the (left) and W (right) pr distributions. 
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FIG. 6. Generated product of the cosine of the angle between the Higgs momentum direction (in the iH centre-of-mass) with 
respect to the iH direction (in the tiH system), and the cosine of the angle of the ^^(left) and (right) momentum (in the 
Higgs centre-of-mass system) with respect to the Higgs direction (in the tH system). 
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FIG. 7. Same as in figure. but after applying the acceptance cuts. 
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FIG. 8. Same as in figure. but using reconstructed objects with truth match. 
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FIG. 9. Same as in figure. but without truth match. 




FIG. 10. Reconstructed product (without truth match) of the cosine of the angle between the Higgs momentum direction (in 
the tH centre-of-mass) with respect to the tH direction (in the ttH system), and the cosine of the angle of the VF^(left) and 
IF”(right) momentum (in the Higgs centre-of-mass system) with respect to the Higgs direction (in the tH system). 




FIG. 11. Reconstructed product (without truth match) of the cosine of the angle between the Higgs momentum direction (in the 
tH centre-of-mass) with respect to the tH direction (in the tiH system), and the cosine of the angle of the Higgs 6-quark(left) 
and t-quark(right) momentum (in the Higgs centre-of-mass system) with respect to the Higgs direction (in the iH system). 
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(Gen.) x^=cos(e“).cos(e|^J 





(Gen.) x^=cos(9'").005(6") 


FIG. 12. Two binned generated product of the cosine of the angle between the Higgs momentum direction (in the iH centre-of- 
mass) with respect to the tH direction (in the tiH system), and the cosine of the angle of the ^^(left) and (right) momentum 
(in the Higgs centre-of-mass system) with respect to the Higgs direction (in the tH system). 
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FIG. 13. Two binned reconstructed product (without truth match) of the cosine of the angle between the Higgs momentum 
direction (in the tH centre-of-mass) with respect to the tH direction (in the ttH system), and the cosine of the angle of the 
^^(left) and £~ (right) momentum (in the Higgs centre-of-mass system) with respect to the Higgs direction (in the tH system). 
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FIG. 14. Higgs reconstructed mass (without truth match) fit (see text for details). The dashed line correspond to the signal 
and shaded region the combinatorial background. 







